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Medium-induced parton energy loss, resulting from gluon exchanges between the QCD matter 
and partonic projectiles, is expected to underly the strong suppression of jets and high-pr hadron 
spectra observed in ultra-relativistic heavy ion collisions. Here, we present the first color-differential 
calculation of parton energy loss. We find that color exchange between medium and projectile en- 
hances the invariant mass of energetic color singlet clusters in the parton shower by a parametrically 
large factor proportional to the square root of the projectile energy. This effect is seen in more than 
half of the most energetic color-singlet fragments of medium-modified parton branchings. Applying 
a standard cluster hadronization model, we find that it leads to a characteristic additional softening 
of hadronic spectra. A fair description of the nuclear modification factor measured at the LHC may 
then be obtained for relatively low momentum transfers from the medium. 



High transverse momentum partons (pr 2> 10 GeV) 
produced in heavy ion collisions interact with the QCD 
matter in the collision region while branching. This in- 
teraction is thought to cause the strong medium modi- 
fication of hadronic spectra and jets measured in heavy 
ion collisions at the LHC and at RHIC. The modeling 
of this jet quenching phenomenon has focussed so far on 
medium-induced parton energy loss prior to hadroniza- 
tion [1[ , assuming that for sufficiently high px hadroniza- 
tion occurs time-delayed outside the medium. However, 
if the color flow of a parton shower is modified within the 
medium, then hadronization can be affected irrespective 
of when it occurs. Here, we analyze for the first time 
the medium-induced color flow in a standard parton en- 
ergy loss calculation. Compared to the current modeling 
of parton energy loss 0-0] , this will be seen to result in 
a characteristic softening of the ensuing hadronization. 
It may thus affect significantly the extraction of medium 
properties from the measured nuclear modification factor 
at the LHC (l^. 

We start by considering an elementary building block 
of a parton shower, the q —¥ q g parton splitting. For 
a small light-cone energy k + of the gluon compared to 
the parent parton, x = k + /p + <C 1, and for transverse 
gluon momentum k with Kn = k/k 2 , the gluon spectrum 
reads, to leading order in a s , 

+ dl vac a s C R v 2 m 

In the presence of QCD matter, interactions between pro- 
jectile parton and medium result in a modified spectrum 
d/ rad = d/ vac +G?/ mod that can be written as an expansion 
in powers of the ratio £ = L + /Ajj of in-medium path 
length L + and an elastic mean free path A^j (for details, 
see [lj]). To first order in this opacity expansion, the 
medium-induced radiation spectrum reads 
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FIG. 1: N = 1 opacity diagrams for gluon radiation from 
a projectile quark in the large- N c limit. The most energetic 
color-singlet clusters are denoted by thicker lines and cor- 
respond to color flow between projectile components (upper 
diagrams, contribution dJJJ^) or between projectile and tar- 
get (lower diagrams, contribution dI{^^ M ). Diagrams on the 
right hand side include a 3-gluon vertex. 

Here, (...) denotes averaging over the transverse mo- 
mentum transfer q of a single interaction, and Ki = 
(k — q)/(k — q) 2 . Medium-induced interference effects 
enter via the factor 

/ s%L+)\ f 1 for 1/ui « L+ , 
1 \ ui^L+ J \0 for 

(3) 

Thus, the medium-modification d/ mod can occur only 
for quanta of sufficiently small formation time = 
2k + 1 (k — q) 2 <C L + . To first order in opacity and at 
large N c , we identify three contributions with distinct 
color flow 
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Here, d/™° t d arises from probability-conserving virtual 
corrections that do change neither color flow nor kine- 
matic distributions in the projectile. The contributions 
dI^° d M and dI^^ M contain medium effects and can be 
characterized by the invariant mass of their leading clus- 
ter, as we discuss now. 



2 



For a projectile with (light-cone) energy much larger 
than that of its target scattering partner, p+ 3> t+ , and 
for x <C 1, the most energetic parton in the final state 
in Fig. [T] is defined unambiguously and carries energy 
p+ = (1 — x)p+ . Distributing the energy of the gluon 
equally between its q and g-legs in the large AT c -limit (our 
final conclusions will not depend on this assumption), the 
energy and invariant mass of the most energetic singlet 
cluster in a vacuum splitting read 

k 2 /2a;. 



(l + x/2)p+ 



(5) 



These relations hold for the leading color singlet clusters 
of the vacuum contribution (|T|), for the virtual correc- 
tion dl™ d in flU, and for dI^° d M (in the latter case, 
the distributions in x and k contain information about 
the medium modification). In contrast, for dI^ d M color 
flows from the most energetic final state parton directly 
to a target component of momentum t, and 
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Here, Qt = V%Eth and Eth is the typical (thermal) en- 
ergy of the target component. The invariant mass in ([6]) 
is parametrically larger than in the vacuum case ([5]) by 
a factor \fp+. This motivates the choice of subscripts in 

jrmcd „ J Jrmed 
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One can determine the kinematic conditions and time 
scales required for dJ^^ M to contribute. The color- 
inclusive sum (j4]) depends on the formation time l/io^ . 
To first order in opacity, we find that the color-differential 
contributions on the right hand side of Q carry inter- 
ference factors involving two formation times and 
1/ujq = 2k+ /k 2 (see ll| for full result). Here, we focus 
on the relevant limiting cases. 

If parton energy loss becomes negligible (w{~ L+ <C 1) 
and if the formation time of the final state gluon is large 
(u)q L+ < 1), then 
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In this case, the color exchange between the medium and 
the projectile occurs predominantly at early times before 
the splitting and hence the color flow and invariant mass 
of the leading color singlet is vacuum-like, dI^ d M = 0. 
In contrast, in the limit 1/cjq <C L+ <C l/^i, one 
finds that dl^ t d ~ -3(K 2 ), dI£* d M ~ (K 2 ), and 
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even if the color-inclusive 



2(Kr/). Thus 
sum ((4]) vanishes, there is a kinematic regime in which 
medium-induced color flow changes the invariant mass of 
the leading color singlet significantly. 

If there is parton energy loss, one finds in the limit 
l/wf , 1/wq < L+ 
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In summary, whenever the color-averaged contribution 
dU to radiative energy loss is significant, the resulting 
medium modification of color flow changes the kinematic 
properties of the most energetic color singlet in more than 
50% of all cases, dl^M > 0.5dl mcd . 



For a gluonic projectile, one finds to first order in 
opacity, and in the large N c limit a similar fraction of 
clusters C tar g of parametrically large invariant mass ([5]) . 
To higher orders in opacity, multiple interactions between 
projectile and target enhance the fraction of color singlet 
clusters Ct arg [HJ- Therefore, the calculation above illus- 
trates generic features of medium-modified color flow in 
models of parton energy loss. 

Since color is conserved in QCD and since hadrons 
are color singlets, a dynamically consistent hadroniza- 
tion prescription must relate hadrons to color singlet 
fragments of the parton shower. We therefore expect 
on general kinematic grounds that the medium-modified 
distribution flSJ) of leading color singlets has consequences 
for the dynamics of hadronization. We now illustrate this 
point in a cluster hadronization model that encodes es- 
sential features of the prescription implemented in the 
MC event generator HERWIG 12j; the parton shower 
is evolved perturbatively to a hadronic scale (Qo ~ 600 
MeV in Ref. [12J) at which the entire shower is decom- 
posed into its color singlet components ('clusters'). Clus- 
ters C of low invariant mass (Aic < -M cr = 4GeV) are 
then decayed directly into pairs of hadrons, while clus- 
ters of larger invariant mass are decayed first into pairs 
of daughter clusters, C -> IV, till the invariant mass 
of the daughters satisfies the condition for decay into 
hadrons. In parton showers evolved in the vacuum, far 
more than 90% of all clusters are found to have a small 
invariant mass, M.c < 4 GeV 12j. We therefore take in 
the following Mc vlLC < 4 GeV. We want to adapt this 
cluster hadronization model to a medium-modified par- 
ton shower, for which at the final scale Qo the most en- 
ergetic parton carries energy p+ . In more than 50% of 
these showers, color will flow from this parton directly to 
the target, and one finds for the most energetic cluster 
Ctarg an invariant mass -Me targ > 4 GeV within the range 
Pf > Per = -Mct/Qt- Estimating the thermal energies of 
target scatterers by their ideal gas value E t h ~ 2.7 T, we 
find for T = 200 (500) MeV a value of Q T ~ 760 (1900) 
MeV corresponding to p+ = 21.0(8.4) GeV. Therefore, 
for a wide, phenomcnologically relevant range of parton 
energies p~j > p+, hadronization of Ctarg involves an ad- 
ditional step Ctarg —> XY that is absent in the fragmen- 
tation of cluster C va c, and that leads to a most energetic 
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daughter cluster X of energy [lj 
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Therefore, amongst the distribution of final clusters C/, 
defined as those that decay directly into hadrons, the 
daughters of C targ are significantly softer (C/ = X, P£ < 
p + f ) than C vac (Cf = C vac , P^ >Pf), see ©. 

To illustrate how the softening of the cluster distri- 
bution due to color flow may affect transverse momen- 
tum spectra, we consider a sample of hadronic collisions 
with final state parton showers evolved to scale Qo. We 
parametrize the resulting transverse partonic spectrum 
by a power-law (within the geometry of a heavy ion col- 
lision, the large momentum component of a hard process 
points transverse to the beam direction) 
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A fraction f t of the partons in (|10p arc endpoints of clus- 
ters Ctarg with color flowing directly to the target. The 
distribution of final clusters Cf resulting from C targ reads 
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Here, A = Qo\J Pa/Qr, and the last factor in (fT2")) 
comes from transforming (fTX)|) from to p£ with the 
help of ((9]). The first term in ([T2"T) comes from clus- 
ters Ctarg of small energy P + < P+ and small invari- 
ant mass, that decay directly into hadrons, Cf = Ctarg- 
The second term comes from clusters Ctarg of large in- 
variant mass that undergo a cluster decay, C/ = X. 
Both mechanisms contribute in an intermediate range 
P+ — A < P + < P+ where one finds an enhancement 
F(P + ) > 1. A more realistic treatment, including correc- 
tions to the limit p~^ 3> t + , may be expected to smoothen 
the O-functions in (fT2j) . Here, we steer clear of these 
mo del- dependent uncertainties by focussing on the high- 
energy region P + > P+ where an additional cluster de- 
cay results in a suppression, F(P + ) < 1. For phenomeno- 
logically motivated parameter choices, namely a power- 
law spectrum (fTU)) with n = 6, values Qo = 600 MeV 
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FIG. 2: The color-induced suppression factor F(P + ) of (fl2|) 
plotted for sufficiently large cluster energies P + > P cr where 
additional color-flow induced cluster decays occur. 



and .Mcr = 4 GeV typical of hadronization models, and 
a range of expected thermal energies, the color-induced 
suppression persists up to cluster energies well above 100 
GeV, see Fig. H 

Medium-effects in single inclusive hadron spectra are 
typically expressed in terms of the nuclear modification 
factor Paa(pt) = (dN AA /dp T ) / (ra coU dN^ / dp T ) ■ Let 
us denote by P^A the nuclear modification factor cal- 
culated in a standard implementation of parton energy 
loss in which the quenching dynamics is factorized from 
hadronization and the latter is treated as in the vac- 
uum [l[. In this case, the energy degradation due to 
d4| is accounted for in P^A , whereas effects of medium- 
induced color flow are not. The latter effects may be 
included by subjecting a fraction f t of the hadronic yield 
to a further suppression (fT2")) at the cluster level. We 
do this by relating the energy of the final clusters C / to 
the hadronic transverse momentum, P + = v2— , and 
approximating the momentum fraction z carried by the 
hadron by the average of a boosted, isotropic two-body 
decay, (z) ~ 3/4, 
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The resulting Paa is plotted in Fig. [3] starting from a 
standard implementation of parton energy loss for central 
Pb-Pb collisions at the LHC (curve for Paa taken from 
Fig. 5 of Ref. [Hj]). We chose this baseline since - in con- 
trast with many results in the recent literature - it refers 
to a relatively small density of the medium (character- 
ized by the quenching parameter q = 1 GeV 2 / fm) , and 
therefore significantly overpredicts the Paa measured at 
the LHC. Supplementing this model with a color-flow 
induced suppression F in f t = 50% (100%) of all frag- 
mentation patterns according to (1131) . we find that Paa 
is reduced significantly by ~ 0.1 (~ 0.2) compared to 
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Baseline R AA '» qhat = 1 GeV'/fm 

Baseline & color flow, t '= 50%, Q T = 760 MeV 

Baseline & color flow, f = 100%, Q T = 760 MeV 

9f ALICE (preliminary) 

© CMS (preliminary) 




FIG. 3: The nuclear modification factor Raa(pt). The base- 
line calculation of kinematic effects (solid black curve) is sup- 
plemented with the effect of color- flow modified hadronization 
according to (fT3")l . 



R l j$. For /( = 100%, this brings the model into reason- 
able agreement [16[ with preliminary results of ALICE 
and CMS at the LHC. 

Recent measurements of the nuclear modification fac- 
tor Raa in central Pb-Pb collisions at the LHC have es- 
tablished unambiguously that Raa(pt) increases notably 
at intermediate transverse momentum px € [5; 20] GeV, 
and they provided first information about the continuing 
slow rise of Raa(pt) in the range pt € [20; 100] GeV. 
This has been used in several recent studies 0-0] to con- 
strain medium properties entering the modeling of parton 
energy loss. It is a common feature of all these models 
to assume gluon exchanges between the projectile par- 
ton and the medium, while neglecting in their dynamical 
implementation the ensuing changes in the color struc- 
ture of the parton shower. Although heuristic models 
of medium-modified hadronization have been discussed 



previously [14l |15| , we assessed here for the first time the 
color-differential information that is implicitly contained 
in standard implementations of parton energy loss. From 
Figs. [2] and [3] we learn that medium-modified color flow 
may affect single inclusive hadron spectra significantly 
up to the highest transverse momenta (pr > 100 GeV, 
say). Also, remarkably, in contrast with the temperature 
dependence of the baseline R^a > the smaller the thermal 
energy in (|12j). the larger the color-flow induced contri- 
bution to the suppression of Raa- As a consequence, a 
significantly smaller density of the medium may then be 
sufficient to account for the observed suppression of Raa ■ 
Medium-modified color flow is an inevitable conse- 
quence of models of parton energy loss. However, its 
numerical manifestation in Raa may depend sensitively 
on the microscopic implementation of parton energy loss, 



including e.g. the spatio-temporal embedding of the par- 
ton shower in the QCD matter and the probability that 
the produced partons escape the medium without gluon 
exchange. As illustrated in Fig. |3l the effect of medium- 
induced color flow is potentially large and thus needs to 
be constrained in phenomenological applications. It is 
thus important to include the effect of medium-modified 
color flow in full microscopic simulations of parton en- 
ergy loss, where also refined or different hadronization 
prescriptions can be explored. Such further studies could 
give insight into the existence or absence of changes in the 
hadrochemistry and in the fragmentation pattern of jets 
in heavy ion collision. They could also help to understand 
any possible difference between the Raa of reconstructed 
jets and leading hadrons. 
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